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ABSTRACT 

Using self-consistent cosmological simulations of disc galaxy formation, we analyse the 
1.4 GHz radio flux from high-redshift progenitors of present-day normal spirals within 
the context of present-day and planned next-generation observational facilities. We 
demonstrate that while current radio facilities such as the Very Large Array (VLA) 
are unlikely to trace these progenitors beyond redshifts z > 0.2, future facilities such as 
the Square Kilometer Array (SKA) will readily probe their characteristics to redshifts 
z < 2, and are likely to provide detections beyond z « 3. We also demonstrate that 
the progenitors of present-day cD galaxies can emit in excess of 10 /xJy of flux at 
redshifts z > 1, and may be a non- negligible contributor to the micro- Jansky source 
counts derived from current deep VLA cm- wave surveys. 

Key words: methods: A^-body simulations — galaxies: disc — galaxies: elliptical and 
lenticular, cD — galaxies: formation — galaxies: evolution — radio continuum: galaxies 



1 INTRODUCTION 

Astronomy is on the brink of an instrumental revolution. 
The next generation of observational facilities will afford 
in most cases, in excess of 1—2 orders of magnitude im- 
provement in sensitivity and/or resolution across the full 
electromagnetic spectrum - from the X-ray {X-Ray Evolv- 
ing Universe Spectrometer), to the optical [Extremely Large 
Telescopes) , to the infrared [James Webb Space Telescope), 
to the millimeter [Atacama Large Millimeter Array), and 
into the radio regime [Square Kilometer Array - SKA, here- 
after). The SKA^ is expected to deliver sub-^Jy sensitivity 
at 1.4 GHz - ~50— 100 times that capable with the existing 
premier radio facility, the Very Large Array (VLA) which 
will be upgraded to the Expanded Very Large Array (E-VLA) 
and whose sensitivity will increase by factor five before the 
SKA will be built. 

At 1.4 GHz, radio counts are dominated by active 
galactic nuclei (AGN) at fluxes > a few mjy, while the 
faint er (~10 AtJy) counts are dominated by starburst galax- 
ies ^Windhorst et al. 1985. 1993: Hopki ns ct al. 2000). At 
the anticipated SKA sensitivity limit (~0.1 liij), little is 
known of the sources which may domin ate the rad i o uni - 
verse. Based on an empirical model, Hop kins et al. I (|200(t ) 
predict that normal disc galaxies near z ~ 0.5 — 1 will dom- 
inate the sub-/xJy counts. As the radio continuum flux is 
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known to be roughly proporti onal to the star formation rate 
(SFR) for non- AGN systems dCondon Ill992h . one might an- 
ticipate that the SKA would probe the progenitors to these 
normal disc galaxies directly, should this prediction be valid. 

To test the above hypothesis, we have constructed cos- 
mological simulations of disc galaxy formation which include 
a self-consistent treatment of temporal evolution of the ra- 
dio flux. Our simulations wer e carried out using the gal ac- 
tic chemodynamics code GCD+ iKawata fc Gibson l2003allbl) . 
GCD+ is a three-dimensional tree V-body/smoothed particle 
hydrodynamics (SPH) code which incorporates self-gravity, 
hydrodynamics, radiative cooling, star formation, super- 
novae (SNe) feedback, and metal enrichment. GCD+ takes ac- 
count of the chemical enrichment by both Type II (SNe II) 
and Type la (SNe la) SNe, mass-loss from intermediate mass 
stars, and follows the chemical enrichment history of both 
the stellar and gas components of the system. The radio flux 
is calculated based upon an empirical law in which the radio 
flux is linked with the event rate of SNe II as desc ribed in 
Section ( Brcssa mTilva fc Granato II2002I : ISilva l lT999) . 

The details of GCD+ relevant to the specific radio source 
count simulations described here are presented in Section |2| 
In Secti0n r3.ll we show the expected fiux from progenitors of 
simulated disc galaxies and compare the result with both the 
current observational data and the anticipated SKA sensitiv- 
ity limit. Section [3.21 described what we consider to be plau- 
sible candidates for some fraction of the faint radio sources 
detected in current faint cm-wavelength radio surveys. Our 
results are discussed fully in Sectional 
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Figure 1. Dark matter density map of a portion of the simulation volume of a comoving 20h~^ Mpc diameter sphere (middle panels), 
and predicted J-band image (physical scale) of galaxies Dl (upper panels) and D2 (lower panels), over the redshift range 2=1.45 to z=0. 
The AB-magnitude system has been employed. 



2 METHOD 

2.1 The Galactic Chemodynamical Evolution 
Code: GCD-^ 

Our simulations were carried out using GCD+, our original 
galactic chemody namical evolution code. Details of the code 
are pre sented in iKawatal 1^9^ and iKawata fc Gibson I 

1 2003a|) . The cod e is e ss entially based on TreeSPH 
Hernauist fc Katzl Il989l : iKatz. Weinberg fc Hernq uist 
W9t ) , which combines the tree algorithm iBarnes fc Hut 
19861) for the computation of the gravitational forces 
with the SPH jLucv I Il977l: iGingold fc Monaghan I Il977^ 
approach to numerical hydrodynamics. The dynamics of 
the dark matter and stars is calculated by the A^-body 
scheme, and the gas component is modeled using SPH. It is 
fully Lagrangian, three-dimensional, and highly adaptive in 
space and time owing to individual smoothing lengths and 
individual time steps. Moreover, it includes self-consistently 
almost all the important physical processes in galaxy forma- 
tion, such as self-gravity, hydrodynamics, radiative cooling, 
star formation, SNe feedback and metal enrichment. 

Radiative cooling which depend s on the metallicity of 
the gas (derived with MAPPINGSHL lSutherland fc Dopita I 



Il993ll is taken into account. The cooling rate for a gas with 
solar metallicity is larger than that for gas of primordial 
composition by more than an order of magnitude. Thus, the 
cooling by metals should not bo ignored in immcrical sim- 
ulations of galax y formation (.Kaellander fc Hultman ..199^ : 
iKav et al. iboodl . 

We stress that GCD+ takes account of energy feed- 
back and metal enrichment by both SNe II and SNe la, 
and metal enrichment fr om intermediate mass stars (see 
iKawata fc Gibson 1120033 . for details). The code calculates 
the event rates of SNe H and SNe la, and the yields 
of SNe II, SNe la and intermediate mass stars for each 
star particle at every time step, considering the ISalpeterl 
(1955) initial mass function (mass range of 0.1^60 M^n) 
and metallicity dependent stellar lifetimes iKodama Ill997l : 
iKodama fc Arimoto I IT997i) . We assume that each mas- 
sive star (> 8 Mq) explodes as a Type II supernova. 
The SNe la rates are calculated using the model pro- 
posed by Kob avashi. Tsuiimoto fc Nomoto (2000). The 
yields of SNe II, SNe la and intermed i ate mass stars 
are ta ken froin^^osi^^^^^^^^_^9^^^^^^^^0 
il999l) and Ivaj^der^Hoe^^^GroeiMwegen 1 Jl997ft . respec- 
tively. The simulation follows the evolution of the abun- 
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dances of several chemical elements (^H, *He,^^C, ^''N, ^^O, 
^"Ne, ^"Mg, ^^Si, ^''Fe, and Z, where Z means the total 
metallicity). GCD+ makes it possible for us to follow the 
galaxy formation processes, taking into account both chem- 
ical and dynamical evolution effects self-consistently, and 
analyse the observational properties, such as optical and ra- 
dio emission, from simulation results. 



2.2 Cosmological Simulations 

Using GCD+, we carried out a series of high-resolution simu- 
lations within an adopted standard A-dominated cold dark 
matter (ACDM) cosmology, ^0=0.3, Ao=0.7, ^2b=0.019/l~^ 
h—0.7 and CT8=0.9 ^. We used a multi- resolution technique 
in order to maximise the spatial resolution within the re- 
gions where the disc progenitors form and evolve. Tidal 
forces from neighbouring l arge-scale structure are in cluded 
self-consistently (see also iKawata fc Gibson I l2003lJ) . The 
initial conditi ons for the si mulations were constructed us- 
ing GRAFIC2 llBertschinger 1.20011 AU our simulations use 
isolated boundary conditions. Initially, we perform a low- 
resolution A'^-body simulation of a comoving 20 h-^ Mpc 
diameter sphere which consists of 137376 particles with the 
mean separation of 20 /i~^/64 Mpc. The mass of each parti- 
cle in our low-resolution simulation was 3.63 X 1O^M0, and 
a fixed softening of 9 kpc (in physical scale) was applied. 

At redshift z — 0,we selected a 10 Mpc diameter spher- 
ical region which contains several galaxy-sized dark matter 
halos. We traced the particles which fall into the selected 
region back to the initial conditions ai z — 43.5, and iden- 
tify the volume which consists of those particles. Within this 
arbitrarily shaped volume, we replace the low-resolution par- 
ticles with particles a factor of 64 times less massive. The 
initial density and velocities for the less massive particles 
are self-consistently calculated by GRAFIC2, taking into ac- 
count the density fields of a lower-resolution region. Finally, 
we re-run the simulation of the entire volume (20 Mpc 
sphere), including gas dynamics, radiative cooling, and star 
formation. The gas component is included only within the 
high-resolution region. The surrounding low-resolution re- 
gion contributes to the high-resolution region only through 
gravity. The mass and softening length of individual gas 
(dark matter) particles in the high-resolution region are 
7.33 X 10'' (4.94 X lO'^) Mq and 1.14 (2.15) kpc, respectively. 

Using the friends-of-friends (FOF) technique, we next 
identified six stellar systems which consisted of more than 
2000 star particles. Two of these systems possess disc-like 
morphologies and internal kinematics consistent with be- 
ing supported primarily by rotation (and have > 10000 star 
particles each). Table summarises the basic properties of 
these disc galaxies which we refer to as galaxies Dl and D2 
hereafter. In Table Q the virial mass Mvir is defined as the 
mass within the virial radius r^h , where the latter is defined 
to be the radius of the sphere containing a mean density 
of 178 f2n '^^ times the c ritical value, pcrit = 3Ho/87rG, fol- 
lowing ^ke^N^arr^X__Frcnk (199^. The rotation veloc- 
ity Viot was measured directly from the rotation curves of 
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Figure 2. Predicted optical V-band magnitude and radio 1.4 
GHz luminosities for galaxies Dl (triangle) and D2 (sq uare). 
Crosse s show the observational data for nearby galaxies from lBell I 
i2003t). 



the respective gas discs. Galaxy Dl has a similar rotation 
velocity to that of the Milky Way, with a Ti'-band lumi- 
nosity approximately twic e that of an L* late-ty pe spiral 
{Mk, = -22.98±0.06 mag: lKochanek et al. 1200111 . Because 
the stellar masses associated with our simulated galaxies 
are likely to be overestimated, due to the ubiquitous over- 
cooling problem (White & Frenk 1991), we consider galaxy 
Dl to be a fair representation of a typical disc galaxy. Galaxy 
D2 is approximately a factor of two more massive and lu- 
minous than Dl. Admittedly, a sample of two is not ideal, 
but we emphasise that this is a pilot study for the longer- 
term goal of better characterising the capabilities of the next 
generation of radio facilities. Since progenitors of these two 
galaxies have similar radio fiuxes despite their different for- 
mation histories, as shown below, it is re-assuring at least 
that our conclusions are independent of the specifics of the 
two high-resolution disc galaxies described here. 

Fig. shows the evolution of the dark matter distribu- 
tion in a central portion of the simulation volume (middle 
row), and the evolution of the stellar component (J-band) in 
a 100-kpc region centred on galaxies Dl and D2 (upper and 
lower rows, respectively). Here, the optical luminosities are 
derived by our population synthesis code adopting the sim- 
ple stellar population of the public spectrum and chemical 
evolution code, KA97 iKodama Ill997l : iKodama fc Arimoto I 



^ These cosmological parameters are close to the most 
favoured va lues after the Wilkinson Microwave Anisotropy Probe 
JSpergel et al. II2OO3I) . The differences are no larger than 15 %. 



I1997D . taking into account the age and me tallicity of star 
particles (see also iKawata fc Gibson ll2003bl and references 
therein). We have taken into account fc-corrections, but have 
neglected the influence of dust extinction. Both galaxies 
form through conventional hierarchical clustering by redshift 
z—1; the disc component is built-up via smooth accretion of 
gas thereafter. By the present-day, both simulated galaxies 
have prominent stellar and gaseous discs. 
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Table 1. Simulated disc galaxy properties. 



Name 

Dl 
D2 



Mvir 

(Mq) 
9.8 X 10" 
2.4 X 10^2 



(kpc) 
305 
408 



Vrot 

(km s-i) 
220 
230 



Mk 
(mag) 
-23.9 
-24.3 



Number of Particles 
Gas dark matter Star 
3997 17408 11819 

12222 41792 27636 



2.3 Analysis: Radio Continuum 

We assume here that a simulated galaxy's radio continuum 
consists of both thermal (L^) and non-thermal {L^"^) emis- 
sion components. We do not currently include a contrib ution 
from AGN. FoUowing lBressan. Silva fc Granalol (120021 ) . the 
non-thermal emission is assumed to follow a simple empirical 
formula related directly to the SNe II rate i?sNo ii, 



L [v) = 



_ RsVie II 

yr" 



+ 1.38(- 

11/ -lTir-l\ 

— (erg s Hz ), 



49GHz 



(1) 



where v is rest-frame frequency (in GHz). Since GCD+ takes 
into account stellar lifetime, i^sNoii is different from SFR. 
We assume that the thermal emission follows the sim- 
ple scaling relation described by ISilva I l|l999h (see also 
ICondon fc YinlllQQOl) . 



10.9 



(ghz) 



(2) 



While admittedly a simplistic assumption, these empir- 
ical laws do reproduce successfully the radio flux of observed 
present-day disc galaxies. Fig.|21shows the predicted optical 
V'-band and 1.4 GHz fluxes for galaxies Dl (triangle) and D2 
(square). The crosses c orrespond to the observational data 
for nearby galaxies from lBelll l|200^. As we do not yet have 
a self-consistent treatment of dust extinction within GCD+, 
we have de-reddened the observational data in Fig. fusing 
the extinction tabulated byt Bell (2003 ) . We assume that the 
luminosity within an aperture of 50 kpc (in physical scale) 
radius encompasses the entire optical and radio-emitting re- 
gions. Fig.|5|shows that the predicted optical and radio lumi- 
nosities for the simulated galaxies are consistent with those 
of observed galaxies. Encouraged by this, below we assume 
that this empirical law is also valid for high-redshift galax- 
ies. In the next section, we discuss the temporal evolution 
of the predicted radio flux for the progenitors of our two 
simulated disc galaxies. 



3 RESULTS 

3.1 Radio flux from disc progenitors 

Fig. |3| shows the predicted 1.4 GHz flux (observer's frame) 
from progenitors of galaxies Dl and D2 at different red- 
shifts. Here, we measure the total flux within an aperture of 
50 kpc (in physical scale) radius and take into account the 
fc-correction. At each redshift, we use the FOF to find the 
stellar systems identified as the progenitors of galaxies Dl 
and D2. Only progenitors which consist of more than 2000 
star particles are displayed in Fig.|3 Kt z = 1.8, galaxy D2 
has two progenitors, which by redshift z = 1.4 have merged. 
Actually, even by visual inspection, galaxy D2 appears to 
have more than two progenitors still identiflable at z — 1.45, 
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Figure 3. The predicted 1.4 GHz flux in the observed frame 
from progenitors of the simulated galaxies Dl (triangles connected 
by doted line) and D2 (squares connected by solid line) at dif- 
ferent redshifts. Cross es correspond to the observational data in 
iHaarsma et al. 1 i200(]l) . The open circle represents the predicted 
flux from a progenitor of a present-day (simulated) cD galaxy in 
our model. Thick dashed line indicates the anticipated SKA flux 
sensitivity limit (0.1 /.iJy) in 12 hours. 



all of which could be reso lved with the SKA arcsecond spa- 
tial resolution'^ (see also ICollev et al'~lll997^ . However, as 
our simple FOF classified them as one system, we simply 
define such system as a single object in the accompanying 
figure. 

The aperture radius of J?ap =50 kpc has been chosen 
arbitrarily. This is roughly a radius which encompasses the 
entire stellar system identified as a single stellar system by 
our FOF at any redshift which we focus on. To see the aper- 
ture effect, we also measured the fulxes within i?ap = 100 
and 10 kpc. At z < 1, the flux within _Rap ~ 100 kpc agrees 
with that within J?ap ~ 50 kpc within about 10 %, while 
the flux within J?ap = 10 kpc is systematically (about 15 % 
in the worst case) smaller than that within i?ap ~ 50 kpc, 
because the stellar systems are bigger than 10 kpc (see Fig. 
0. At 2; > 1, since the central galaxy is relatively small, the 
contribution from the satellite galaxies are not negligible. 
Consequently, _Rap ~ 100 kpc leads to 50 % larger flux than 
J?ap = 50 kpc in the worst case. On the other hand, in the 
multiple merger like the progenitors of galaxy D2 at z = 1.45 
(Fig. 0, i?ap = 10 kpc sometimes leads to about 20 % of 
flux within i?ap ~ 50 kpc for the biggest galaxy. Note that 
the radio interferometer, such as the SKA, enables to re- 
solve them, and they are likely detected as multiple sources, 
as mentioned above. 



10 kpc at redshift z = 
the adopted cosmology. 



1.45 corresponds to ~1.2 arcseconds in 
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Our quoted SNe II rates are averages over the most 
recent 1 Gyr for any system at z < 1 and 0.05 Gyr for 
the ones ai z > 1. While a somewhat arbitrary "binning", 
we confirm that our results are not sensitive to the specific 
bin-size adopted. Such a sampling suppresses the stochastic 
appearance of relatively rare short period (~Myr) starbursts 
which are likely caused by resolution limits. Thus, the results 
of Fig. 13 represent a mean trend in the expected flux from 
progenitors of galaxies Dl and D2, as implicitly supported 
by the smooth nature of the curves in the accompanying 
figure. 

The c rosses in Fig. El corresp ond to the observed VLA 
dataset of iHaarsma et ahl 11200(1) . The flux limit for this 
dataset is approximately 10 /iJy. Our simulations suggest 
that the VLA count source data corresponds to the (recent 
- z < 0.2) progenitors of present-day disc galaxies. This is 
consistent with the conc lusions drawn by optical fo l low u p 
for the same sources by iRoche. Lowenthal fc Kool l)2002il . 
At higher-redshifts, these progenitors possess radio fluxes 
significantly below the feasible 1.4 GHz fiux limit of the 
VLA. Their predicted fiux reaches ~1 ^Jy at redshift z ~ 1, 
and decreases to ~0.1 /iJy by redshift 2 = 2 (with little evo- 
lution at redshifts z > 2). We do not trace the progenitors 
beyond redshift z > 3, as the number of bound star particles 
become too few (<2000) for the results to remain reliable. 
Galaxies Dl and D2 show very similar trends for the tem- 
poral evolution of the observed 1.4 GHz fiux. Based upon 
our (admittedly limited) sample of two high-resolution sim- 
ulated disc galaxies, we conclude that the trends shown by 
the curves in Fig. |3 are fairly representative of typical disc 
systems. Consequently, our results suggest that the SKA is 
capable of detecting the progenitors of normal disc galaxies 
at redshifts z <3, provided it meets its design specifications 
of 0.1 /iJy (the expected 5 a lev el achievable after a 12 hour 
integration, e.g. iHopkins et al. 1 12000) as shown by the thick 
dashed line in Fig.|3 

Fig. 21 shows the SFR of the progenitors of galaxies Dl 
and D2 at different redshifts. Here, the quoted SFR has been 
derived by averaging over the same aperture and temporal 
bin as that described previously. The cross es represent the 
inferred SFR of the VLA faint sources from lHaarsma et al. I 
l|200(]() . who estimates the SFR using a n empirical rel ation 
between SFR and radio continuum fiux dCondon 1199211 . Be- 
cause of the adopted empirical link between radio fiux and 
SNe II rate (eq. 0), the latter of which also relates directly 
to the SFR, Figs. 13 and 0] provide somewhat complementary 
information. The reason why there are no observational data 
with low SFR at high redshift is the observational selection 
effect due to the limited sensitivity, which is the same as the 
effect which leads to the well-known envelopes of observa- 
tional data in the absolute radio luminosity-redshift plane 
(e.g. Fig. 9 of lGilbert et al. Il2004^ ■ Fig. gl suggests that the 
progenitors of galaxies Dl and D2 have much smaller SFRs 
than the VLA faint sources at z > 0.2 (which is simply an- 
other way of stating why they have much lower radio fiuxes) . 

The dotted and solid histograms of Fig.^Jshow the his- 
tory of the total SFR for galaxies Dl and D2, respectively. 
Here, the total SFR means the summation of the SFR of 
star particles which fall into the region within a galactocen- 
tric radius of 50 kpc at redshift 2 = 0. For both galaxies, the 
total SFR is similar to the SFR of the most massive progen- 
itor at redshift z < 1.5. However, at z > 1.5 the SFR of the 



o 
o 





O 






























o 


X X 


x^ 





o 


r X >kx 






; \ X 








: y,>^ X 








01 










redshift 



Figure 4. The SFR for progenitors of the simulated galaxies Dl 
(triangles connected by doted line) and D2 (squares connected 
by solid line) at differ ent redshifts. Crosses correspond to the 
observational data from lHaarsma et al. Dotted and solid 

histograms show the history of the total SFR for galaxies Dl 
and D2 (see text for details). The open circle shows the SFR for 
a progenitor to a present-day (simulated) cD galaxy. The reason 
why there are no observational data with low SFR at high redshift 
is the observational selection effect due to the limited sensitivity. 



most massive progenitor is significantly lower than the total 
SFR. This difference indicates that a significant number of 
stars are forming in smaller progenitors which later accrete 
onto the most massive progenitor. An obvious conclusion 
from this cursory analysis is that in order to properly esti- 
mate the observed fiux from high-redshift objects, it is (not 
surprisingly) important to take into account the hierarchi- 
cal clustering history of the systems. Galaxy D2 is a clear 
example - at z > 1.5, we have noted previously the presence 
of two massive progenitors, and the fact that we are able to 
analyse the radio flux from each progenitor separately. This 
demonstrates that numerical simulations can be a powerful 
tool for studying the anticipated radio flux from individual 
galactic building blocks. 



3.2 Radio flux from cD galaxy progenitors 

In the previous section, we found that the VLA sources ob- 
served at z > 0.2 are unlikely to be progenitors of normal 
disc galaxies. The obvious question then arises as to which 
population they belong. We suggest that one candidate 
might be the progenitors to bright central-dominant (cD) 
cluster galaxies. Using the same multi-resolution technique 
described in Section 12.21 we carried out a high-resolution 
simulation of a cluster of galaxies which was similarly iden- 
tified first from a low-resolution cosmological simulation of 
a comoving 40 Mpc diameter sphere, and whose virial 
mass at redshift z = was 1.2 X 10" M0. The mass and 
softening length of individual gas (dark matter) particles in 
the high-resolution region were 5.86 x 10' (3.95 x 10*) Mq 
and 2.27 (4.29) kpc, respectively. 

The J-band image of the simulated cluster redshift 
z = 1.3 is shown in Fig. |3 The basic properties of the clus- 
ter at this redshift are summarised in Table |5| The cluster 



6 D. Kawata, B.K. Gibson & R.A. Windhorst 



-100 





X (kpc 



100 



Figure 5. J-band (in tlie observed frame, and using the AB- 
magnitude system) image of a progenitor of a cluster of galaxies. 



has a massive central galaxy - a progenitor of what will be- 
come a present-day cD galaxy - and a number of associated 
member galaxies. The 1.4 GHz flux and SFR within the cen- 
tre 50 kpc (in physical scale) of the cluster are denoted as 
open circles in Figs. |3| and |3 As seen in Fig. |S| the flux 
within this region is dominated by the flux from the progen- 
itor of the cD galax;y itself. Figs. Eland 0] demonstrate that 
this cD galaxy progenitor has a sufBciently high SFR, and 
consequently a sufficiently high 1.4 GHz flux, to reach the 
present-day VLA detection limit. Therefore, we suggest that 
the progenitors of present-day cD galaxies may be lurking 
within the faint VLA source counts. Fig. |3 shows that the 
SFR of this cD progenitor is ~300 Mq yr'^. Since SFRs of 
this magnitude are likely to be accompanied by associated 
dust, the intrinsic ultraviolet and optical light would suffer 
from consequent extinction, re-appearing as sub-millimeter 
emission. Thus, the VLA faint sources which are also de- 
tected b y the Sub-millimeter C ommon User Bolometric Ar- 
ray (e.g. lChapman et al. lEool) may be counterparts of our 
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axy. ims result is consistent witn tne pre- 
. 'Governato et al. Ill998l: iNagamine Il2002t 
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2004) which demonstrate that high red- 



shift starburst galaxies, such as Lyman break galaxies, likely 
evolve into large galaxies in clusters and groups at z=0. 



4 DISCUSSION AND CONCLUSIONS 

Combining our cosmological simulations with an empirical 
law for radio continuum flux emission has enabled us to trace 
the temporal and spatial evolution of radio flux from star- 
forming disc and central dominant cluster galaxies. We have 
found that current radio telescopes, including the VLA, are 
capable of tracing the 1.4 GHz flux of the progenitors of 
today's disc galaxies only locally [z < 0.2). Progenitors of 
normal disc galaxies have a 1.4 GHz flux <10 /iJy ai z > 0.2, 
with their predicted flux reaching ~0.1 /iJy at z = 2 — 3. 
Since the flux limit of the SKA at 1.4 GHz is expected to 
be ~0.1 fiJy in 12 hours, it should be capable of detecting 



disc galaxy progenitors at redshifts z < 3. In addition, from 
a cluster simulation, we also found that progenitors of cD 
galaxies can possess a 1.4 GHz flux in excess of 10 fiJy at 
2 > 1, and may be a candidate for some fraction of the faint 
VLA sources at 2 > 1. 

It is worth reminding the reader that hierarchical clus- 
tering often leads to the presence of more than one signifi- 
cant progenitor (or building block) for any one disc galaxy 
at high redshift. Our simulated galaxy D2 provides a typ- 
ical example of such a system. We predict that the radio 
flux from these individual building blocks will be observ- 
able by the SKA. At the sensitivity limit probed by the 
SKA, these building blocks of normal galaxies may dominate 
the radio sour ce counts, as they d o in the optical faint blue 
galaxy counts (jPriver et al. Il99^ . Observations of such ex- 
tremely faint radio continuum sources may therefore provide 
a unique window into the accretion history of galaxies. In 
addition, the SKA may be able to provide redshift i nforma- 
tion for the sources through HI line observations (e.g. lBraun I 
1996). Then, the SKA would be capable of providing the lu- 
minosity function of radio sources at different redshifts (e.g. 
ICowie. et al. 2004), which would give strong constraints on 
disc galaxy formation scenarios. To interpr et such observa- 
tions, the oretical models are im portant (e.g. lTakeuchi et al. I 
[2OOI; Na gashima et al~ll2002^ : our pilot study shows that 
self-consistent cosmological simulations will provide useful 
theoretical support for such analyses. Encouraged by our 
preliminary results, we are undertaking a series of larger 
volume simulations (anticipated completion in 2005), in or- 
der to improve the statistics associated with the program 
described here. 

One note of caution is that the current numerical 
simulatio ns suffer from the we ll-documented over-cooling 
problem l| White fc Frenki I199J) . which results in an un- 
avoidably high SFR in the high-redshift building blocks, 
and a consequent underestimation of the SFR at later 
epochs. In tandem, this leads to disc systems which suf- 
fer from having overly massive stellar halos and bulges, 
at the expense of a prominent thin disc. Improved mod- 
eling of SNe feedback is believed to be crucial in over- 
coming this long-standing problem, and it has been sug- 
gested that stronger SNe feedback than that adopted here 
may lead to end-products w hich more closely resemble ob- 
served disc galaxies today (Thacker & Couchman 200 ll 
Igonim cr-Larscn. G 5 tz_fc Portinari 2003; Robertson ct alj 
I20Q4I iBrook et al.1f2004h . Having said that, even models 
which do incorporate strong feedback only result in SFR 
histories a factor of two diff erent from the cano nical models 
presented here (e.g. Fig. 2 of'Broo k et al. l2004^ . Hence, this 
particular problem unlikely impact our predicted radio flux 
by more than a factor of two. Nevertheless, more realistic 
numerical models of galaxy formation remain the constant 
target for our group (and others), and are a necessity for 
contributing to the design and anticipated performance of 
next-generation facilities such as the SKA. 



ACKNOWLEDGEMENTS 

We are grateful to Nobuo Arimoto and Tadayuki Kodama 
for kindly providing the tables of their SSPs data. We ac- 
knowledge the Yukawa Institute Computer Facility, the As- 



Cosmological simulations of the high-z radio universe 7 



Table 2. Properties of a cluster progenitor at z = 1.3. 

Mvir Tvir a"' M Number of Particles 

(Mq) (kpc) (km s^^) (mag) Gas dark matter Star 
6.3 X 10" 544 185 -28.2 59365 140947 70793 
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